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Mass transfer in eccentric binaries: the new Oil-on- Water SPH 
technique 
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ABSTRACT 

To measure the onset of mass transfer in eccentric binaries we have developed a two-phase 
SPH technique. Mass transfer is important in the evolution of close binaries, and a key issue is 
to determine the separation at which mass transfer begins. The circular case is well understood 
and can be treated through the use of the Roche formalism. To treat the eccentric case we use 
a newly-developed two phase system. The body of the donor star is made up from high- 
mass water particles, whilst the atmosphere is modelled with low-mass oil particles. Both 
sets of particles take part fully in SPH interactions. To test the technique we model circular 
mass-transfer binaries containing a 0.6 M© donor star and a 1 white dwarf; such binaries 
are thought to form cataclysmic variable (CV) systems. We find that we can reproduce a 
reasonable CV mass-transfer rate, and that our extended atmosphere gives a separation that is 
too large by aproximately 16%, although its pressure scale height is considerably exaggerated. 
We use the technique to measure the semi-major axis required for the onset of mass transfer 
in binaries with a mass ratio of ^ = 0.6 and a range of eccentricities. Comparing to the 
value obtained by considering the instantaneous Roche lobe at pericentre we find that the 
radius of the star required for mass transfer to begin decreases systematically with increasing 
eccentricity. 

Key words: gravitation - hydrodynamics - methods: numerical - binaries: close - stars: 
evolution - stars: mass-loss - X-ray: binaries 



1 INTRODUCTION 



Mass transfer in binary stars is a very important process that 
completely changes their evolution. Several eccentric binaries that 
are undergoing mass transf er exist. One such system is the low- 
mass X-ray binary Cir X-1 jMurdin et alj|l980l : IXauris etalJI 19991 : 
I Johnston et~Zl ll999': 'Clarkson et al."2004l). Its orbital parameters 
were analysed by Tauris et al. ( 1999) who argue that the system is 
probably a 2 Mq star orbiting a neutron star with an eccentricity of 
about 0.9. Higher mass systems are more frequently observed; one 
such example is the Be/ X-ray transient A0538-6 6, located in the 
Large Magellanic Cloud jKretschmar et al.ll2004 l). This consists of 
a B2 Ille star orbiting a neutron star with a period of 16.65 d in 
an eccentric orbit with e ^ 0.7. In most of these systems the mass 
transfer between the stars is driven by winds rather than through 
Roche lobe overflow. 

Systems such as Cir X-1 are thought to be created when the 
heavier star explodes in a supernova. Another possible way to cre- 
ate a hard eccentric binary is through a close encounter between a 
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circular binary and a single star or binary. The rate of such events 
is only significant when the stellar density is high, so they occur 
primarily in environm ents such as the cores of globular clusters 
dOavies & Benz|[r995l) . For a density of about 10' stars/pc^ the av- 
erage time between close encounters for a system is about I Gyr. 
Hence most of the binaries in a dense globular cluster core undergo 
at least one encounter. 

The theory of mass transfer in the circular case, where the 
flow is continuous and steady, has been studied for a lone time and 
is well understood jPaczvnsH197ll : lRenvoize et alj2002h . For two 
stars in a circular orbit there exists a corotating frame in which 
they are stationary and if the stars are taken to be point masses the 
potential, known as the Roche potential, is well defined. The first 
connected equipotential surface that surrounds both stars is known 
as the Roche lobe and provides a boundary for the potential well in 
which the star sits. If we make the assumption that the envelope of 
the star deforms to the equipotential surface whilst conserving its 
volume then we can define the Roche lobe radius such that a sphere 
with that radius has the same volume as the Roche lobe (Eggletoi] 
Il983h . This can then be used in conjunction with a spherically sym- 
metric model of the stellar structure, such as that produced by a 
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stellar evolution code, to predict the point of onset of mass transfer. 
The eccentric case is much more difficult because the assumptions 
made in order to derive the Roche potential are no longer valid. The 
relative velocities and distances of stars in an elliptic orbit vary so 
there is no corotating frame. 

A commonly used method to assess whether mass transfer oc- 
curs in an eccentric binary is to make the assumption that the stars 
corotate at periastron and apply the circular theory there. How- 
ever it is not clear that it is possible to define an instantaneous 
Roche lobe because, among other things, to do so assumes that the 
timescale upo n which the star adjust s to the ch anging force is suf- 
ficiently short ICharles et al. 1983; B rown & B oyle 1984). 

Another approac h is to simulate such systems numerically. 
lBovle&Walke3 ( [T986) ma de simulations with t est particles placed 
on the surface of the donor. lHavnesetal]h98d) simulated Cir X- 1 
with a core and an extended atmosphere of test particles. How- 
ever they only looked at a single orbit. Full numeric al simula- 
tions of complete systems are rare but lRegosetai] ( l2005l) presented 
some low-resolution smoothed particle hydrodynamic (SPH) sim- 
ulations where they considered four different eccentricities for a 
single semi-major axis a and mass ratio q. 

In this paper we introduce a new technique (called Oil-on- 
Water) within the SPH formalism. We define two types of SPH 
particles, heavy water particles that make up the stellar interior and 
very light oil particles that sit on top of the star. We are thus able 
to resolve the mass transfer even though the fraction of the stellar 
mass transferred per orbit is very small. We apply this technique 
to mass transfer in eccentric binaries and investigate a variety of 
binary systems, varying the semi-major axis a and eccentricity e. 

In Section 2 we discuss eccentric mass-transferring binaries 
from a theoretical point of view, as well as the observational evi- 
dence for their existence. In Section 3 the oil-on-water method is 
described and some tests presented. The results of the simulations 
are given in Section 4 and discussed in Section 5. 



2 THE OIL-ON-WATER TECHNIQUE 

In order to investigate mass transfer in eccentric binaries a two- 
phase SPH technique has been developed. This section describes 
the new approach and some tests of it. 



2.1 Smoothed particle hydrodynamics 

SPH was developed by ILucvI ( Il977h and iGingold & MonaghanI 
( Il977h . It has been widely used in various astrophysical applica- 
tions. SPH is a particle-based Lagrangian scheme: the particles' 
motions follow the fluid velocity. This differs from grid based meth- 
ods where the fluid flow between a grid of cells is measured. An ad- 
vantage of SPH is that effort is not spent solving the hydrodynamic 
equations in regions of space that are devoid of matter; in the case 
of a binary star system simulated in a grid-based code this is most 
of the computational volume. Another feature of SPH is that one 
can follow t he ev olution of the part icles. For reviews of SPH see 
iBenj ( Il990h and iMonaghanl ( Il992l) . Here only a few key points 
will be mentioned in order to show how our implementation differs 
from the standard technique. Our code is based upon the code of 
[Benz ( 1990). 

At the heart of SPH lies the kernel W(r, K) and the smoothing 
length h. These set the size and shape of a particle's sphere of influ- 
ence; how much a particle affects other particles that lie at position 
r. For example the density at r is given by 



p(r) = £m,lV(|r-r,|,/i). 



(1) 



One simply takes the sum of all the particles' masses weighted by 
W{\r — rj\, h). If particle j is far from r its contribution to the sum 
is negligible (or even zero), whilst if it is close it has a large contri- 
bution. 

SPH is a Lagrangian averaging scheme so it is easy to take 
derivatives and hence simple to calculate the force on a particle /. 
For example the derivative of the density is given by 



{Vp(r))= V ^VW(k-r,|,ft). 
^ Pi 



(2) 



Use of a similar expression for the derivative of the pressure P and 
symetrisation leads to the acceleration of a particle /, 



dv, 

— = — > nij 



p, p-j 



VW(r,/i)- VO, 



(3) 



where Pi is the pressure at particle i; the first term is the sym- 
metrized expression for the pressure gradient. The artificial viscos- 
ity Yijj is introduced in SPH to improve th e treatment of shocks. 
We us e the standard formulation for this bv lMonaghan & Gingold 
( Il983l) . We utilise a simple polytropic equation of state with y = 
5/3. Finally, VO gives the gravitational force owing to all the other 
particles. For a star in stable hydrostatic equilibrium the gravita- 
tional force on a particle is balanced by the pressure gradient. 

Gravity is not a local force so all particles have to be taken 
into account and hence a brute-force calculation scales as A'^ in 
the particle number. This would mean that the method would lose 
much of its usefulness because a large amount of computational 
effort would be used for the gravitational force summation. We use 
instead the hierarchical tree method of iBenz et al. I h990l) which 
employs a binary tree and scales as A' log A'. 

To deal better with the large range in densities that arise 
in sim ulations of these kind we employ a variable smoothing 
length teenz|[l990l) . The derivative of the smoothing length is cal- 
culated as 



At ^ 



(4) 



We modify the derivative thus calculated to keep the number of 
neighbours that each particle has between 80 and 120, subject to a 
maximum A of 0.1 R^. 



2.2 The two phase technique: oil-on-water 

A typical mas s-transfer rate fo r a low-mass close binary is 10"^ to 
lO"' Mg yr"' ( |Pattersoniri984l) . This implies that, in a typical orbit 
of a few days to a fraction of a day, about 10"'" to 10"'^ Mq is 
transferred. To resolve this mass flow using equal mass particles 
we would require at least 10'^ of them and this is beyond what can 
be accomplished today. Within SPH it is possible to use a range of 
particle masses but if the mass range is large numerical problems 
arise. As an example we can consider the density calculation of 
Equation [T] If a particle of unit mass comes just within the range of 
a light particle with a mass of lO"'' the heavy particle dominates the 
density completely, even when the kernel at that point is very small. 
We circumvent this problem by introducing a two-phase scheme. 
The interior of the star is made up of heavy water particles while 
the atmosphere contains very light oil particles. It is these particles 
that take part in mass transfer. We also introduce an artificial force 
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Figure 1. A cartoon figure of the force balance between the number gradient 
force f „ and the gravitational force Fg,av The number gradient of the heavy 
water particles introduces a force perpendicular to the surface of the star In 
order for this force to be smooth a sufiicient number of water particles is 
needed. This is accompHshed, in pail, by a remapping of the water particles, 
see Section lT4l 



to keep the two different types of particles apart. This is calculated 
from the number gradient rather than from VP. 

In our simulations we use 15,390 water particles to simulate a 
0.6 Mo star so the average mass of a water particle is 3.92x 10"' Mq. 
The mass of each oil particle is 10"'* Mq so the average water par- 
ticle has the same mass as 4 x 10' oil particles. We use 32,691 oil 
particles in our runs, so the total mass in oil particles is negligible. 

Figure[T]illustrates how the oil particles are balanced on top of 
the water star. We have introduced an artificial force, Fp, to sepa- 
rate the oil layer from the water particles. This is based on the num- 
ber density of the water particles and therefore is perpendicular to 
the edge of the water surface. This prevents the oil particles from 
penetrating into the stellar interior and balances their gravitational 
attraction to the star. It is defined according to 



F„ = CkV«„, V«, = V V, W(|r - r^l, KJ 



i 



(5) 



where the sum is over the water particles j. To control the force we 
introduce two parameters, a and p. They are defined as 



^GMgrn,- 



(6) 
(7) 



To calculate V«, the smoothing length of the water, /t„ is used rather 
than the usual average between the particles. A small value of a 
means that the oil and water particles must come closer to one an- 
other before they begin to feel the repulsive force and so the oil 
layer lies closer to the star. On the other hand the force must change 
smoothly otherwise a typical oil particle will come towards the wa- 
ter particles with a large velocity, bounce on the hard force barrier 
and get a large velocity kick away from the star. A particle's time 
step decreases when the rate of change of the force is large and this 
slows down the code. The strength of the force is controlled by /? 
and the same arguments apply. Optimal values for the parameters 




Figure 2. The left graph shows a schematic depiction of the number density 
of water particles computed using the SPH method. The derivative of this is 
shown to the right. If an oil particle comes within |/! of the water edge the 
repulsive force from the water particle density gradient declines and it can 
fall in towards the centre of the stai'. 



Force on 
water particle 



Force on 
oil particle 



Force from water particle Standard 
Force from oil particle Standard 



Additional force F„ 
Standard 



Table 1. The different force laws applied when calculating the force on a 
particle due to one of the others. Standai'd denotes the normal SPH force 
calculation. 



have been deduced from a large number of empirical tests. For the 
production runs presented in this paper we use 



and 



P = 40. 



The force balance on the particles is then 



dv, 

dr = -z.' 



Pi p] 



VW(r,/2) - — - VO, 

nii 



(8) 



(9) 



(10) 



where the number gradient force is only applied for interactions 
between oil and a water particles: see Table[T]for further clarifica- 
tion. The Fn force is similar to the normal pressure force (see Equa- 
tion[3]l except that the acceleration it causes has no dependence on 
the mass of the particles. For interactions between two oil particles 
or two water particles the ordinary equations are used. This is very 
important because the oil particles are not test particles. They obey 
the full SPH formalism so we can, for example, study the formation 
of an accretion disc around the accreting star. 

2.3 The oil-water kernel 

At the boundary between the oil and water particles the number 
density of the water particles drops quickly (see Figure [2] left 
panel). This gives rise to a problem with the number gradient force. 
Consider two oil particles located at positions A and B in the right- 
hand panel. If the particle at position B is perturbed slightly in- 
wards the number density force increases and the particle is re- 
pulsed. However, if the particle at A moves inwards the number 
gradient decreases and the infall continues. Therefore once a parti- 
cle is within |/i of the stellar surface it falls into the star. This can 
be prevented by keeping the oil particles sufficiently far away from 
the water edge. However this creates a separation between the oil 
particles and the star. 

Another approach to this problem is to change the kernel for 
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the oil-water interaction when calculating F^. The kernel used to 
calculate the number gradient force does not have to be the same 
as that used for the rest of the calculations. Since this interaction 
is very different the demands on this kernel are different. The most 
commonly used kernel is 




^ V < 1 

1 < V < 2 
otherwise 



(11) 



where v = rjh jMonaghan & Lattanzic|[l985l) . Use of this kernel 
causes particles to clump together at a separation of r = ^h. This is 
because the derivative of the kernel VW(r, K) has a minimum at this 
point. In our case this would imply that the oil layer is from 
the edge of the water star. Because we aim to have the oil layer as 
close to the star as possible we investigated three other kernels, 

• a kernel wit hout a break point (m inimum in its derivative), 

• the kernel of iHeranteTd] ( Il994f ) and 

• a kernel Wp(r, h), constructed with a pair of polynomial func- 
tions, with a break point at vj, = 0.4. 

The kernel without a break point proved unsatisfactory be- 
cause as particles come close to one another the timestep falls 
sharply. To prevent this a large value of c\ must be used to prevent 
the oil particles from approaching the star. Hence the oil particles 
are placed far from the water surface. The other two kernels were 
both well-behaved but Wp(r, h) seemed somewhat more robust and 
hence was used for our production runs presented in Section|4] It is 
given by 



W„{r,h) = — 




< V < 1 

1 < V < 2 . (12) 
otherwise 



To obtain this kernel, we first choose it to be the same as Wml 
for r > h, subject to a normalisation factor. We construct a new 
fourth-order polynomial for r < h, subject to the constraints that 

• the kernel and its first and second derivatives are continuous 
at V = 1, 

• the first derivative of the kernel is zero at v = 0, 

• the first derivative of the kernel contains only one extremum, 
the break point, which is a minimum, and 

• the break point is located at Vbreak < f (we choose Vbreak = 0.4). 

These conditions specify the kernel W,, up to the standard normali- 
sation factor. 



2.4 Remapping of the water particles 

It is advantageous to have as many water particles as possible in 
the outer parts of the star and hence close to the interface of the 
oil and water particles. This can be accomplished by remapping 
the water particles to increase the number density at the edge of 
the star. The star is initially built up using a close-packed lattice 
with the water particles placed equidistantly. Particle masses are 
assigned t o obtain the correct den sity profile, following the YREC- 
models o f iGuenther et"ai] l ll992h . To increase the concentration of 
water particles at the surface we remap their positions according to 
a power law of index y = 0.9. Each particle / is updated according 
to 

''/.remap ^norm'", (13) 



where c„orn, is a normalisation factor to keep the stellar radius con- 
stant. This increases the number density at the outer edge by a fac- 
tor of about 1.25 with respect to the centre. The larger number of 
water particles at the edge makes the surface smoother; this allows 
us to place the oil layer closer to the surface. 



2.5 Production of the single-star models 

To build a star the water particles are first set up using the procedure 
described above. The resulting model is then relaxed for approxi- 
mately 100 Ticiax. where T,.eiax is the dynamical timescale given by 



^relax 



(14) 



Then the oil layers are added on top of the relaxed water star. We 
have tried different numbers of oil particles: one and three layers 
of oil require 9,983 and 32,691 oil particles respectively. These are 
added according to a close-packed lattice with a spacing equal to 
the internal smoothing length, h = 0.025 Rq. The oil particles each 
have a mass of 10"'* Mq. The model is relaxed until the particles 
have settled into a stable configuration; this typically takes about 

50 Trelax- 

The two stars with one and three layers of oil are shown in 
Figure [3] In the model that originally had a single layer the oil par- 
ticles are evenly spread around the star but there is no depth in the 
oil layer. The other has visible structure. The oil particles all have 
the same mass so density is well-mapped by the number density 
of the oil particles. There are more oil particles close to the water- 
oil boundary so the atmosphere is densest here. The simulations 
presented here all utilise 32,691 oil particles as in the three-layer 
model. The remapping of the water particles is also visible in the 
density contrast between the centre and surface of the water particle 
star. 



3 CIRCULAR BINARIES 

In order to test the reliability of the oil-on-water model we used the 
code to model circular binary systems. Because the potential field is 
circular, co-rotating binaries is well described by the Roche formal- 
ism there is a large body of work that describes them and hence use- 
ful comparisons may be made. Two questions can be posed. First, 
can we make a binary system that reproduces observed properties 
correctly? Second, how much does the unphysically large size of 
the atmosphere in our models affect the results that we obtain? 
We answer these two questions by making models of a cataclysmic 
variable system in the manner described below. 



3.1 The Roche potential 

In a circular binary where one star is filling its Roche lobe, the mass 
flow is continuous and steady, and a number of simplifications can 
be made. By assuming that the stars are centrally condensed for 
the purposes of calculating the potential and that they are corotat- 
ing w ith the orbit we obtain the Roche potential jPringle & Wadd 
Il985l) . In a corotating frame centred on star 1, with the line con- 
necting the two stellar centres along the x axis, the potential takes 
the form 



GM, 



^{x - a)2 + / + ^2 
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X (Rq) X (Rg) 

Figure 3. Two examples of a relaxed composite star. The particles plotted are those within h of the xy plane. The small dots are water particles and the large 
dots oil particles. In both models there are 15,390 water particles. In the left-hand model there are 9,983 oil pailicles: in the right-hand 32,691. In the right-hand 
figure a clear density structure can be seen in the oil particles. 



(15) 



where a is the separation, SIq the angular velocity and the reduced 
mass yu = M2l(M\ + Mi). The last term in this expression is due 
to the centrifugal force. The stationary point between the two stars 
is the inner Lagrangian or LI point, and the surface around each 
star that passes through this point is its Roche lobe. The radius of a 
spherical star with the same volume as its Roche lobe is called the 
Roche-lobe radius i?L and is fitted within a few percent by 



a 



QA9q 



2li 



0.6^2/3 -Hln(lH- ^1/3)' ^^^^ 

( lEggleton|[T983l) . where q - M2/M1 is the mass ratio. If one of 
the stars fi l ls its Roche lobe mass flows through LI to the other 
star. lRittejfTgS^ showed that the amount of mass transferred in a 
binary system depends on the pressure scale height Hp of the donor 
star's atmosphere and the difi^erence between the donor star's radius 
R and R]^, according to 



(17) 



where Mq is the mass-transfer rate of a binary that just fills its 
Roche lobe. For a low-mass main-sequence star Mq =^ 10"** Mq yr"' 



and Hp lO"* Rq. 



3.2 Production of circular binary models 

The transition from the spherical potential produced by the self- 
gravity of a single star and the Roche potential owing to an orbiting 
binary system disrupts the structure of our relaxed star unless care 
is taken. For this reason it is necessary to relax a binary system once 
the star has been placed in its orbit before allowing mass transfer 
to proceed. To do this we first add a point mass to the star to rep- 
resent a compact companion. In all these runs we have used a star 
with mass 0.6 Mq and a 1 Mq point mass. Initially we place the star 
and point mass in a frame rotating about the centre of mass with 



a damping force in place to allow the star to adjust to the gravita- 
tional field of the point mass. The star is taken to be co-rotating: it 
is stationary in the rotating frame other than for oscillations during 
relaxation. Once the star has relaxed the model and point mass are 
transformed back into an inertial frame. The star and the point mass 
are both given the correct velocities to place them on the desired or- 
bit. Relaxation takes approximately lOTreiax- 

For an eccentric mass-transfer binary the relaxation takes 
place in a circular orbit at the apocentre separation where no mass 
transfer takes place. In contrast for a circular orbit the separation 
of the stellar centre and point mass is constant and hence there is 
no point in the orbit at which the relaxation can take place without 
mass transfer occurring. To avoid this problem we relax the star in a 
circular orbit that is slightly wider than the widest orbit in which ap- 
preciable mass transfer takes place. After relaxation the separation 
of the orbit is reduced by 0.5% each dynamical timescale, changing 
the positions and velocities of the particles. This provides a series 
of circular models at different separations that can be evolved fur- 
ther to model a series of binaries. 

At low mass-transfer rates we encounter numerical problems. 
Because the rate of flow of particles into the accretor's Roche lobe 
is small there are only a few particles present, which causes the 
pressure forces to vary very rapidly. These particles gain large ve- 
locities and are ejected from the Roche lobe. To counter this prob- 
lem we turn off sph forces for oil particles with fewer than 50 neigh- 
bours: that is, such particles move balistically. 



3.3 Results 

All the circular binaries consisted of a 0.6 Mq star, modelled using 
the two-phase oil-on-water technique and a 1 Mq point mass. Such 
binaries are believed to be relatively high-mass cataclysmic vari- 
able (CV) systems. In these hydrogen-rich material accretes on to 
a white dwarf and bums explosively at the surface. Ongoing mass 
transfer is thought to be driven by angular momentum loss owing 
to either a magnetic wind or gravitational radiation. Observations 
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Figure 4. The mass-transfer rate in our circular binary models. Tlie crosses 
show the rate at which mass is accreted in the circular' systems and the 
dashed line is a straight-line fit to these poi nts. The so l id line is the observed 
accretion rate for CVs from Equation 17 o f lPattersonI fl984l) . and the dotted 
lines its 1-sigma error interval. 



of many such systems have been made and their accretion rates can 
be measured from their luminosities. We ran models at a range of 
separations, which yielded a corresponding range of mass-transfer 
rates. We measure the rate by fitting the number of particles in the 
Roche lobe of the accretor as a linear function of time, using a stan- 
dard least-squares fitting algorithm. The accretion rate for a given 
binary was roughly constant until a substantial fraction of the enve- 
lope had been transferred. The accretion rates we obtain are shown 
in Figure O along with an observa tional measurement of the CV 
accretion rate from Pattersonl jl984h . 

Mass transfer from a CV is driven by a loss of angular momen- 
tum from the system and the mass-transfer rate is seen to be a func- 
tion of the binary period. At the points where the observed accretion 
rates and the accretion rates obtained from our code are equal we 
can take the separation to be that predicted for the CV we are mod- 
elling. We find this equilibrium separation to be a = 2.03^[5 jjj Rq. 
Applying Equation [16] with q = 0.6 gives a Roche lobe radius 
R]^ = 0.68^q55P Rq. The mean volume of the star, as defined by 
the volume occupied by the water particles, was estimated by find- 
ing its extreme points along the three natural axes of the system 
and calculating the volume of an ellipsoid with those axis lengths. 
The radius of the sphere with the same volume, in analogy to the 
definition of the Roche lobe radius, is 0.57 Rq. If we assume that 
the real edge of the star coincides with the outer surface of the wa- 
ter particles this gives an error in the effective radius of the stellar 
model owing to the extended atmosphere of 15 to 18%. 

The mass-transfer rate in our models varies in an approx- 
imately exponential fashion with separation. Hence the model 
shows the qualitative behaviour predicted bv .RitteJ jl98S) . which 
suggests that the atmosphere is sufficiently resolved. By rearrang- 
ing Equation [T7] and substituting for the Roche lobe radius from 
Equation [16] with ^ = 0.6 we obtain 



, M\ R a 
^'^^\tt] = W -03356- 



(18) 



\Mo/ "p 

Equating the final term to the variation in accretion rate from Fig 



Table 2. A table of the runs with different semi-major axes a and eccen- 
tricities e. In all cases the mass ratio q = 0.6. Column 4 indicates whether 
no mass transfer (N), mass transfer (Y) or massive mass transfer (M) was 
observed. For the mass-transfer case the number of particles escaping the 
donor per orbit is given in Column 5. 



a/Re 


e 


''peri /Rq 


Mass transfer 


Oil particles lost 
from donor per orbit 


2.5 


0.0 


2.500 


N 


- 


2.5 


0.15 


2.125 


Y 


30 


2 5 


20 


2 000 


Y 


200 


2 5 


25 


1 875 


Y 


500 


3.0 


0.0 


3.000 


N 


- 


3.0 


0.15 
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ure|4] gives an effective pressure scale height of 0.022 Rq, unsur- 
prisingly very similar to the smoothing length h. The pressure scale 
height of a low-mass main-sequence star is about 10-"* Rq, roughly 
two orders of magnitude smaller than in our models. This means 
that the onset of mass transfer in our models will be much gentler 
than that seen in real systems. 



4 ECCENTRIC BINAIUES 

To simulate eccentric binaries we carried out a grid of runs with dif- 
ferent a and e. In each case the donor was a main-sequence star of 
M = 0.6 Mq and the accreting companion was a white dwarf with 
M = 1 Mq. To relax the star the point mass was initially placed at 
a distance aapo from the star in a circular orbit as explained above. 
When the system had relaxed sufficiently we released it from co- 
rotation and introduced eccentricity in such a way that it was re- 
laxed at apocentre: the spin of the star was set to coincide with 
the angular velocity of the system there. Therefore the star initially 
co-rotates at apastron. The spin rate is not found to change signif- 
icantly during a complete run. Figure|5]is an example of a system 
with a = 2.5 Rq and e = 0.2 which undergoes mass transfer only at 
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Figure 5. An example of mass transfer in an eccentric binary with a = 2.5 Rq and e = 0.2. Here for clarity only oil particles within h of the xy-plane are 
plotted. The position of the point-mass companion is marked with a cross. The mass of the star is O.6M0 and that of the point mass M = 1.0 M©, as in all the 
simulations in this paper. In this binary mass is transferred only near periastron, which occurs when the star is to the left of the point-mass along the ,t-axis. 
This image sequence starts just before periastron passage and, as the star approaches periastron the mass transfer begins; however there is a noticeable lag. 
The mass transfer actually begins just after the closest approach and continues for almost 1/4 of the orbit. 



periastron, so the mass transfer turns on and off in each orbit. One 
such mass-transfer event is shown. 

Table|2]enumerates parameters for the grid of runs completed, 
stating whether mass transfer takes place and, if so, the number 
of oil particles escaping the donor star per orbit. More eccentric 
systems have a narrower range of a in which stable mass transfer 
takes place. This is expected because the periastron separation is 
given by Op = a(l - e), so a small change in e leads to a large 
change in the periastron separation. 

An example of the mass transfer in one of our runs is given in 
Figure [6] The periodic behaviour of the mass transfer is apparent. 
On each orbit mass begins to flow from the donor at pericentre. The 
number of oil particles lost from the donor at each periastron pas- 
sage is roughly constant after the first few orbits. This indicates that 
the mass transfer is stable. One periastron passage of this system is 
also pictured in Figure |5] 

Our models are plotted in the (a, e)-space in Figure |7] Open 
circles denote systems without mass transfer, filled circles stable 
mass transfer and crosses denote that the there is a massive over- 
flow of particles with water particles being transferred. The predic- 
tion of Equation [16] for an instantaneous Roche lobe at pericentre 
is also plotted. Such an approach assumes that the timescale of the 
mass transfer is sufficiently fast to be able to adjust to the instanta- 
neous distance and rotational velocity. The typical orbital passage 



12000 




Figure 6. The cumulative number of oil particles released from the donor 
against time in the case of an eccentric system (a = 2.5 Rq, e = 0.2). Af- 
ter the first initialisation, where a damping force is active, there is a stable 
flow of particles at periastron and this gives raise to the stepped function. 
The lower curve shows the number of particles captured by the companion. 
Most of the particles transfeired are not captured but escape through L2. As 
discussed in section[TT]the potential at L2 is almost the same height as the 
inner Lagrangian point Ll for e = 0.25 
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Figure 7. Results from runs at different a and e. Open circles denote binaries 
where no mass transfer was seen, filled circles are runs with stable mass 
transfer, and crosses mark models where a massive overflow of particles 
was observed and water particles began to transfer. The curves show the 
separations derived for the onset of mass transfer under the assumption that 
the star exactly fills its Roche lobe at periastron. The three curves are for 
stellar radii of 0.6, 0.7 and 0.8 Rq, with 0.8 Rq being the topmost curve. 
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Figure 8. The ratio rp(e)/rp(e = 0) as a function of the eccentiicity. The 
separation required for onset of mass transfer for diff'erent eccentricities is 
rp(e) The line is a linear fit to the points, given by /(e) = 1 - 0.16e. 

shown in Figure|5]shows that this is not the case. The mass transfer 
does not begin until just after periastron passage and continues for 
almost 1/4 of an orbit thereafter. This is due to the time taken for 
the material to respond to the changing potential as the star goes 
round the periastron passage and the time taken for matter to flow 
from the donor to the companion. The instantaneous Roche lobe 
construction also assumes corotation, which only can be true at one 
point in the orbit. 

In order to quantify the effect of the eccentricity we express 
the periastron distance required for onset of mass transfer, rp(e), as 

rp(e) = rp(e = 0) X /(e), (19) 

where rp(e = 0) is the separation required for the circular case 
and /(e) is a correction function to be determined. If the instan- 
taneous Roche lobe formalism was applicable, i.e. the onset of 
mass transfer was determined by the Roche lobe at periastron, we 
would obtain /(e) = 1. Based on our runs in Table [2] we have es- 
timated the eccentricity for which the mass transfer begins for a 
given semi-major axis a. The results are plotted in Figure [8] We 
have extrapolated the data to e = to obtain an estimate of the pe- 
riastron distance at which mass transfer starts in the circular case 
as rp(e = 0) = 2.212 R© . The required distance for onset of mass 



transfer between the stars can be seen to decrease with increas- 
ing eccentricity. The line is a linear fit to the points which yields 
/(e) = 1 - 0.16e. Equation|I9|expressed in terms of the semi major 
axis a gives 

I -0.16e 

amt(e) = amt(e = 0) X — , (20) 

1 - e 

where flmt(^) is the required semi major axis for mass transfer and 
"ratCe = 0) can be obtained from Equation[T6] Note however that in 
our simulations all runs are corotating at apastron and ^ = 0.6, so 
the formula given above is only tested for this case. 



5 CONCLUSIONS AND FUTURE PROSPECTS 

The simulation of mass transfer in eccentric binaries is non-trivial 
because the mass transferred in each orbit is a tiny fraction of the 
total mass. In this article we have presented a technique that treats 
this problem by introducing a two-phase SPH formalism. Very light 
particles make up the outer part of the donor star, while the inner 
part is formed from heavier particles. Thanks to this oil-on-water 
model have we been able to simulate mass transfer in eccentric bi- 
naries. The simulations presented here cover a large range of semi- 
major axes and eccentricities at a mass ratio of 0.6. 

We have found that the onset of mass transfer in our simu- 
lations does not follow the prescription of the circular case. We 
have measured the eccentricity required for mass transfer given a 
semi-major axis a. Using this we show that the minimum distance 
between the stars that leads to mass transfer decreases linearly with 
eccentricity. 

The oil-on-water model could be used to simulate other in- 
teresting astrophysical processes. One example is the onset of 
common envelope evolution. A close binary system undergoes a 
common envelope (CE) phase if the secondary star is unable to 
accrete all the mass transferred by the donor jPaczvnskH Il976l : 
lRasio&Liviolll996h . The core of the evolved star and the sec- 
ondary orbit one another in a cloud of gas, which is eventually ex- 
pelled. This transports angular momentum away from the system 
and the period of the system decreases. The oil-on-water technique 
can follow the angular momentum transferred at the onset of the 
common envelope phase accurately, making it an useful tool for 
such an investigation. 
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